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ABSTRACT 



Comet 17P/Holmes underwent the largest cometary outburst in recorded his- 
tory on UT 2007 Oct. 23, releasing massive quantities of dust and gas. We used 
the Canada- Prance-Hawaii telescope to obtain wide-field images of 17P/Holmes 
on 15 dates over a period of 3 months following the outburst and employ them 
here to examine the subsequent activity of the nucleus and the nature of the 
ejecta closest to the nucleus. Through aperture photometry we observed the in- 
ner coma (within 2500 km of the nucleus) to fade from an apparent magnitude 
of 11.7 mag to 17.6 mag, corresponding to absolute magnitudes of 8.1 mag and 
12.4 mag, between UT 2007 Nov. 6 and 2008 Feb. 12. A second much smaller 
outburst occurred on UT 2007 Nov. 12, three weeks after the original outburst, 
suggesting that the nucleus remained unstable. The surface brightness profile 
of the inner coma was consistently shallow relative to the expected steady-state 
profile, and showed a persistent brightness enhancement within ~ 5000 km of 
the nucleus. We propose that sublimating ice grains created an ice grain halo 
around the nucleus, while fragmenting grains were responsible for the shallow 
surface brightness profile. 

Subject headings: comets: individual (17P/Holmes) 
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Introduction 



Comet 17P/Holmes underwent the largest cometary outburst in recorded history on 
UT 2007 Oct. 23.7 ± 0.2 (iBuzzi et al.||2007|), brightening from ~ 17*'^ magnitude to ~ 



2nd ^2 hours (iMoreno et al. 20081 Li et aL 20111). The Jupiter family comet is known 



to have experienced at least two prior outbursts: the first led to its discovery on UT 



1892 Nov. 6 (Holmes 1892) while the second occurred in January 1893 (Denning 1893). 



Whipple (1984) speculated that the unusual outburst activity might have been triggered by 



a grazing encounter with a natural satellite. In this scenario, the second outburst would 
have been caused by the final impact of the satellite with the nucleus. A third major 
outburst practically eliminates the likelihood of this scenario, suggesting that an alternative 
mechanism causes these unusual periods of activity in 17P/Holmes. The two largest 
outbursts - those that occurred in 1892 Nov. and 2007 Oct. - both took place approximately 
five months after 17P/Holmes passed through perihelion at 2.1 AU, suggesting that the 
activity is thermally-triggered, perhaps with a lag due to the time needed for heat to reach 



subsurface ice deposits (Sekanina 2008) 



The mass of material ejected has been estimated to be between 10^° kg and 10^^ kg 



(Sekanina 2008; Reach et al. 2010 Li et al. 2011), corresponding to ~ 0.05% to 5% of 



the nucleus mass. This ejected mass was dominated by the contribution from the dust 
coma that was observed to expand around the nucleus. The peak dust production rate 



approached 10^ kg s"^ on UT 2007 Oct. 24.5 ± 0.01 (Li et al. 2011), about 10^ times 



the dust production rates typical of Jupiter Family comets when at similar heliocentric 
distances. 

In addition to massive quantities of dust, 17P/Holmes released significant quantities 
of volatiles with II2O production rates exceeding 4.5 x 10^^ s~^ (Dello Russo et al. 2008 



Schleicher 2009). Results from large-aperture (radii > 10,000 km) narrowband photometry 
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in the weeks following the outburst led Schleicher (2009) to conclude that enhanced 



water production rates were primarily due to prolonged nucleus activity, with a possible 
contribution from sublimating icy grains in the inner coma. 

In this work, we use a sequence of high-resolution wide-field continuum images to 
examine the evolution of 17P/Holmes during the four months following the outburst. By 
utilizing a single instrument and observing strategy we minimize the effect of systematic 
errors and are able to use a uniform procedure for obtaining, reducing, and analyzing the 
data. This results in a consistent, high-quality data set. In contrast to an earlier study 
of the integrated light from the comet (Li et al. 2011), we here study material near the 
nucleus, which is the most-recently ejected and the most pristine of the observable ejecta. 



2. Observations and Data Reduction 

We used the 3.6 m Canada-France-Hawaii Telescope (CFHT) atop Mauna Kea with 



the wide-field imager MegaCam (Boulade et al. 2003) to obtain high-quality, time-series 



imaging data. MegaCam has a 0.96° x 0.96° field of view provided by an array of 36 
charge-coupled devices (CCDs). Each CCD has 2048 x 4068 pixels, with a pixel scale of 
0".185 pixel^^. We obtained at least one set of five dithered images in the SDSS r' filter 
(Ac = 6250 A, bandwidth = 1210 A) on each of 15 nights between UT 2007 Nov. 6 and UT 
2008 Feb. 12 UT. Throughout 2007 Nov. and on UT 2007 Dec. 11 we obtained two sets of 
exposures with differing exposure times to allow us to probe the bright central region while 
still being sensitive to the faint outer coma. In this work we use the set of exposures with 
the longest exposure time but no saturation of the central region. We have analyzed all 
sets of the exposures and find the agreement between sets obtained on the same night to be 
very good, with differences in calculated magnitudes to be less than 0.02 mag. 
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The observations were taken in queue-scheduled mode. In this mode, non-sidereal 
tracking was unavailable and the comet trailed across each image. Exposure times varied 
between 5 s and 200 s, and were chosen to balance the need for sensitivity while minimizing 
trailing. Table ([T]) provides a journal of observations. 



All images were automatically pre-processed using the Elixir pipeline (Magnier & 



Cuillandre 2004), which uses bias frames and twilight flatfields to remove the instrumental 



signature, and calculates instrumental zero-points using Landolt standard star fields 



(Landolt 1992) observed on each night. We used SWarp, software released by the Terapix 
data center at the Institut d'Astrophysique de Paris, to resample the images and attain 0".2 
astrometric accuracy. The dithered images from each night were median-combined to cover 
the 80" gaps between CCDs and remove cosmic rays (Figure [T]). Each combined image was 
convolved with a two-dimensional Gaussian function with a full-width half-maximum of 
1" to minimize the effect of variations in seeing from night to night. A panel showing the 
central square arcminute of each image is shown in Figure g. The same central regions 
are shown with isophotal contours in Figure ^ to better highlight the evolution of the 
coma morphology. The inner coma initially appears elongated in the anti-solar direction 
but becomes more symmetrical at later dates. 

The sky transparency was monitored using the CFHT Skyprobe and we found that 
only observations made on UT 2008 Feb. 12 experienced significant variations in attenuation 
(> 0.1 mag). For all other nights the weather was photometric. Systematic errors due 
to flatfielding or scattered light are well-characterized for MegaCam and do not exceed 



5% (Magnier & Cuillandre 2004). The sky background was determined in a region as far 
from the observable coma as possible within each stacked image, using the median signal 
within a ~ 4'x 4' box. Although we sought regions on the image devoid of any noticeable 
brightness gradient from which to estimate the background level, it is likely that some coma 
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contamination persisted in our sky background estimations. Specifically, the gas coma was 



observed to expand at a velocity of 1 km s~^ (Drahus et al. 2008), meaning that the gas 
coma reached the projected edges of the MegaCam field of view by 2007 mid-November. 
Fine dust grains dynamically well-coupled to the gas and ejected at high velocities could 
contaminate the outer regions within a month of the outburst, although our observations are 
less sensitive to the smallest grains since their scattering efficiency is small. In Figures (|4]) 
and ([5]) we show the effects of a 10% error on the sky background combined with the 
known sky transparency and systematic errors. Plausible sky subtraction uncertainties 
are negligible in most of our data, given the high surface brightness of the coma. Only in 
observations from January and February 2008 do sky subtraction errors rival systematic 
uncertainies due to atmosphere transparency and calibration uncertainties. 



3. Photometry of the Inner Coma 

The coma was initially observed as an approximately circular dust cloud with a bright 
condensation south-west of the nucleus. As it expanded, the coma became more asymmetric 
with a parabolic nose building in the solar direction and a diffuse tail in the anti-solar 
direction. The nucleus region was always identifiable in the inner coma and stars could be 
seen even through the innermost coma, implying a low optical depth. We used aperture 
photometry to study the material closest to the nucleus with the aim of monitoring the 
comet's continued activity after its massive outburst. 

The position of the nucleus was determined to sub-pixel accuracy using a two- 
dimensional derivative search and used as the center for all aperture photometry. Since the 
geocentric distance increased from 1.62 AU on UT 2007 Nov. 6 to 2.62 AU on UT 2008 
Feb. 12 (Table [T|, an aperture of fixed angular radius would cover an increasing physical 
distance on the sky, and, correspondingly, more of the coma. For this reason, we adjusted 
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the angular radii of the photometry apertures used for each night to compensate for the 
increase in geocentric distance to 17P/Holmes, and instead used apertures of fixed physical 
radius at the comet. Circular apertures of projected radii between 250 km and 25000 km in 
increments of 250 km were used. 

We chose not to subtract the flux of the nucleus from these measurements. The 
expected flux contribution from the nucleus can be calculated using the observing geometry, 



a nucleus radius of 1.62 km, and a geometric albedo of 0.04 (Snodgrass et al. 2006 Lamy 



et al. 2004). The calculated nucleus contribution was typically less than 1% of the flux 
measured within a 2500 km radius aperture. In the worst case, the nucleus contribution 
amounted to 3% of the measured flux but in all cases it was significantly smaller than other 
sources of error. 

Trailing losses are negligible for observations between UT 2007 Nov. 6 and UT Dec. 
11 since the motion of the comet during each exposure was less than the full- width-half- 
maximum of the image. Later images were trailed by up to 7.7 pixels (1".4). To investigate 
the effect of trailing on our photometry we took an image with minimal trailing of 0.2 pixels 
(0".4) and shifted and added duplicate images to mimic trailing. We found that the loss of 
light from a 1.3" radius aperture (2000 km for the image most badly affected by trailing) 
was on the order of 2% and fell to less than 0.1% for a 6". 3 (10^ km) radius aperture. We 
conclude that trailing has a negligible effect on the results from aperture photometry in our 
data. 

We converted fluxes to apparent magnitudes, m, using the instrumental zero-points 
calculated by the Elixir pipeline, correcting for both airmass and the color of the filter. 
To correct for changes in observing geometry we used Equation ([T| to determine absolute 
magnitudes, H - the magnitudes objects would have if observed at heliocentric and 
geocentric distances of 1 AU and a phase angle of 0°. 
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H = m-5 logioirnA) - 2.5 logw^a 



where th and A are the hehocentric and geocentric distances in AU respectively and 
is a phase function that corrects for the dependency of the brightness on the phase angle 
(a). The phase function of the observed dust coma is difficult to measure given the need for 
accurate photometry of the dust over a wide range of phase angles when the coma is bright 
enough to detect. We use an approximation of the phase function for IP/Halley based on a 



fit to data presented in Schleicher et al. (1998) (Li et al. 2011): 



2.5 /o^io$a = 0.045a - 0.0004^2 (2) 

The phase corrections that result from Equation ^ match well with other commonly- 
used phase functions since the change in phase angle was relatively small (< 8°). We 
acknowledge that this is a slight misuse of the concept of absolute magnitude given that 
Equation ([T| refers to unresolved point sources rather than distributed sources such as 
cometary comae. An angular aperture will subtend a different distance at 1 AU than at 
other heliocentric distancs, resulting in an altered amount of material within the aperture. 
However, since we used apertures of fixed physical distances from the nucleus we expect 
this effect to have little impact on our results. Table (|2| shows the results from aperture 
photometry for each night during our observing campaign. 

The surface brightness was calculated from the total flux within concentric annuli 
centered on the nucleus and divided by the area measured in square kilometers. The flux 
per unit area was converted to absolute magnitudes per million square kilometers at the 
distance of the comet, which accounted for changes in the observing geometry during the 
course of our observing campaign. Figure ^ shows the surface brightness proflles after the 
images were convolved with a 2-dimensional 1" Gaussian function. 
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4. Results and Discussion 

4.1. The Fading of Comet 17P/Holmes 

The absolute magnitude within a 2500 km radius aperture increases from 8.1 mag 
on UT 2007 Nov. 6 to 12.4 mag on UT 2008 Feb. 12, indicating a ~ 50-fold decrease in 
brightness (Figure |4]). Both the initial fading observed throughout 2007 November and 
the fading during the entire observing campaign can be reasonably well fitted with an 
exponential function of the form H = a + be'^* where H is absolute magnitude, t is time since 
the first observation in days, and a, b, and c are constants. The best fits, as determined via 
a non-linear least-squares fit, are shown in Figures Q and ([s]). We note that there is no 
clear physical reason for the decrease in brightness to be fitted with an exponential function 



in log-linear space. However, Schleicher (2009) found that similar functions matched the 



changing production rates for OH, CN, NH, C2, C3, and the continuum over a period of 
several months. 



The expansion of the dust coma was studied by numerous observers and the expansion 



velocity was found to be ~ 550 m s ^ (Gaillard et al. 2007 Sekanina 2008 Snodgrass et al. 



2008 Hsieh et al. 2010), although Montalto et al. (2008) reported a lower value of 200 m 



s ^. For this paper we use the relationship between grain velocity, f (m s ^), and radius, 



(/im), derived by Whipple (1951) and assume that a 1 /im radius dust grain has a velocity 
of 550 m s~^ so that: 



(3) 



The high measured velocities imply that the majority of dust particles ejected would 
exit even a 25000 km radius aperture in ~ 1 day. Even 1 cm sized particles (a^ = 10'^) would 
leave a 25000 km radius aperture in about 50 days by Equation ([s]). Given the brightness 
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observed for several months following the outburst, the nucleus must have continued to 
supply the inner coma with fresh material after the 2007 Oct. outburst. This result initially 



seems to be in conflict with the production curve calculated by Li et al. (2011), who found 
that significant activity persisted for only a few days, after which relatively little new 
material was injected into the coma. However, their analysis refers to the entire coma of 
17P/Holmes, which spanned several hundred thousand km, whereas our work focuses on 
a very small fraction of the coma close to the nucleus, and is consequently much more 
sensitive to small, localized changes. 

The apparent magnitude, m, can be used to estimate the total cross-section of the dust 
within an aperture using: 

a = ^-^^ TT rl A' lO-0-4(— ™0) (4) 

where a is the total geometric cross-section of the material (m^), gx is the geometric 
albedo, and rriQ is the apparent magnitude of the Sun. All other variables have the same 
definitions as in previous equations. The variables vh and A are given in AU. The albedo 
of the ejected material was measured to be between 0.03 ± 0.01 and 0.12 ± 0.04 at a phase 



angle of 16° (Ishiguro et al. 2010). The corresponding geometric albedos - as would be 
observed at a phase angle of 0° - are estimated to be 0.05 and 0.20 after using Equation g 
to calculate a phase correction. 

Given the uncertainties, we adopt a nominal value of 0.1 for the geometric albedo of 
dust grains and accept that the actual value may differ from this by a factor of two. The 



apparent magnitude of the Sun in the SDSS r' filter is -26.95 mag (Ivezic et al. 2001). We 
then convert the cross-section into a crude mass, assuming that spherical dust grains with 
radii of 1 /xm and a bulk density of 1000 kg m^ are scattering the light. Figure ([7| shows the 
change in cross-section and mass of the material within a 2500 km radius aperture between 
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UT 2007 Nov. 6 and UT 2008 Feb. 12. The cross-section of material contained within the 
aperture decreases from 6.6 x 10^ m^ to 1.3 x 10® m^. These cross-sections correspond to 
effective circular radii of 46 km and 6 km and masses of 8.9 x 10^ kg and 1.7 x 10^ kg, 



respectively, which was only a tiny fraction of the total mass ejected (2 - 90 x 10^° kg; Li 



et al. 2011). The mass estimates are strongly dependent on the assumed size distribution of 



the particles, and are used here for order of magnitude estimates only. 

For small apertures we can derive mass loss rates using the aperture crossing times of 
the dust grains. Large apertures will contain a mix of old, slow-moving large grains and 
fresh, high-velocity small grains, while the reflected light contained in small apertures is 
dominated by new material that may be only hours old. We continue to use our earlier 
assumptions that micron-sized particles traveling at 550 m s~^ dominate the cross-section 
of observed material. We infer a mass loss rate of 4900 kg s^^ at the beginning of our 
observations, which had decreased to 90 kg s~^ by UT 2008 Feb. 12. Peak mass loss rates 



of 3 X 10^ kg s ^ and 6-7 x 10^ kg s ^ were estimated by Li et al. (2011) and Sekanina 



(2008), respectively. Our results are lower by several orders of magnitude, indicating a 



rapid decrease following the initial outburst. 



4.2. A Second Outburst 

The fading inner coma of 17P/Holmes brightened unexpectedly on UT 2007 Nov. 12 
by ~ 25% (Figure [5]). The expected absolute magnitude, as determined from a fading 
function fit to the November data of the exponential form discussed previously, is 9.3 mag, 
whereas the absolute magnitude calculated from the observed brightness was 9.0 mag. 
Given the associated photometric uncertainties, this is a 3.0 a (99.8% probability) result. 
We note that this significance increases to 3.7 a (99.97% probability) if the UT 2007 Nov. 
12 data point is excluded from the fit. A background object bright enough to cause the 
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increase would need to have an apparent magnitude of 14.6 mag; an inspection of the STScI 
Digitized Sky Survey showed that there was no such object near 17P/Holmes. We attribute 
the brightening to a second, albeit much smaller, outburst on the surface of the nucleus. 
The subsequent observation of 17P/Holmes on 2007 Nov. 13 also lies above the fit shown 
in Figure ([s]), suggesting that some of the material released during the second outburst 
remained in the aperture at this time. 

The surface brightness profiles of UT 2007 Nov. 11, 12 and 13 provide independent 
evidence for a second outburst (Figure [6]). The coma within ~ 5000 km of the nucleus 
is brighter than expected on UT 2007 Nov. 12 when compared to the surface brightness 
profiles of 17P/Holmes on the previous and following nights. We can estimate the onset 
time of the second outburst if we assume that material ejected was dominated by dust 
grains traveling at ~ 550 m s^^. We find that this material could have traveled 5000 km 
in 2.5 hours, leading us to conclude that the outburst started at approximately 2007 Nov. 
12.4. It seems likely that the outburst shut down within hours, given that subsequent 
observations revealed no unusual activity. 

The difference between the expected magnitude and the observed magnitude, as shown 
in Figure ([s]) using a 2500 km radius aperture, amounts to 0.3 mag. Using Equation Q 
and assuming that the material consists of dust grains with radii of 1 /im, albedos of 0.1, 
and densities of 1000 kg m~^, the extra material has a scattering cross-section of 6 x 10^ 
m^, and a mass of 8 x 10^ kg, corresponding to a sphere of material with a radius of 6 m. 

The short duration of the Nov. 12 event suggests that the loss of material did not 
cause significant exposure of volatile-rich patches on the surface of the nucleus. Ejection 
of a surface layer by gas pressure from within the nucleus would likely have left an icy 
crater exposed that would have remained active for days or weeks as volatile deposits 
diminished. The event could have involved the dislodging of a fragment that was only 
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weakly bound to the surface. The fragment would then have broken apart into a collection 
of particles with a much greater total cross-section. Rotational spin-up upon ejection, or 
the sublimation of volatiles that were holding the fragment together could have caused the 
break-up. The nucleus may have remained unstable for a significantly longer period of time 



given that Miles (2009) detected a 0.85 mag brightening on UT 2009 Jan. 4. It is likely 
that 17P/Holmes will continue to undergo minor outbursts for the next few orbital periods 
as unstable regolith settles and freshly-exposed volatiles sublimate. 



4.3. Surface Brightness Profile 

The distribution of material in the inner coma can be studied using the surface 
brightness profile. We fitted the gradient of each surface brightness profile between 10000 
km and 25000 km using a power law. The gradient, rrip, is given by 



d In BiR) 

where B{R) is the flux per 10^ km^, and R (km) is the distance from the nucleus. A 



dust coma generated by an isotropic, steady source would have a slope rup = -1 (Jewitt 



1991). Instead, the measured slope varied between -0.2 and -0.3 with no significant temporal 
trend (Table |3]). We can be confident that the shallow slopes are not a result of inaccurate 
sky subtraction caused by coma-contamination. If the estimated sky value actually included 
a contribution from the dust coma, the images would be over-subtracted, resulting in an 
artificially steep slope. The shallow slope suggests that either the dust production rate 
varied in such a way as to temporarily produce a smaller scattering cross-section close to 
the nucleus, or that the cross-section increased with increasing nucleo-centric distances, 
possibly through a process such as dust fragmentation. The likelihood of the activity on 
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the nucleus varying with time in such a way as to consistently produce a shallow surface 
brightness profile for several months is extremely low. Another process is needed to account 
for the shallow surface brightness gradients. 

Cometary ejecta is likely to include agglomerate grains, which are clusters of smaller 



particles that are loosely held together by electrostatic forces or volatile material (Clark 



et al. 2004). Heating of these tenuous aggregates by incident solar radiation could cause 



them to break apart, greatly increasing their scattering cross-section. For example, the 
brightness of a 1 cm radius cluster of micron-sized dust grains could increase by a factor of 
~ 10^ if it broke apart into its constituent particles. Thus, fragmenting or disaggregating 
particles could cause a shallow surface brightness profile where the scattering cross-section of 



material increases with distance from the nucleus. Yang et al. (2009) identified populations 



of both hot and cold grains in the inner coma of 17P/Holmes using infrared spectra. They 
postulated that the cold ice grains and hot dust grains could have originally been contained 
within a porous aggregate that broke apart when exposed to sunlight. We conjecture that 
fragmenting aggregates were responsible for the shallow profile at small nucleo-centric 



distances (< 25000 km), Hsieh et al. (2010) suggested that this process continued out to 
distances of greater than 10^ km. 



4.4. Ice Grain Halo 

The surface brightness profiles consistently show a bump in the innermost region before 
the profiles straighten out to a decreasing slope (Figure [6]). To investigate the nature of this 
bump we created a model surface brightness profile for each night with a slope matching 
that measured at large angular radii in the CFHT data, plus an unresolved central source. 
The model was forced to match those brightness levels observed at 250 km and 10000 km, 
and then convolved with a 1" Gaussian function to simulate the effects of seeing. Figure (Isl) 
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compares the results of our model (dashed line) to the data (solid line). The comparison 
shows that the brightness of the central region cannot simply be attributed to a smooth 
coma with a contribution from the nucleus. One possible interpretation is that a halo of icy 
sublimating grains persisted around the nucleus. Ice grain halos in comets were suggested 



decades ago (Delsemme & Miller 1970) and evidence for halos has been reported in several 



comets since then, including C/1996 B2 Hyakutake (Harris et al. 1997) and 103P/Hartley 



2 (Meech et al. 2011). Schleicher (2009) found that derived water production rates for 



17P/Holmes during the month following the outburst were too large to be explained by 
sublimation of water ice from the surface of the nucleus. He concluded that a population of 
icy grains must have been sublimating near the nucleus, which fits well with the ice grain 
halo scenario. 

The extent of an ice grain halo is dependent on the sublimation rate of the ice grains, 
which in turn depends on several factors including the size, albedo, temperature, and 
density of the grains. To explore these factors, we calculate the sublimation rate assuming 
thermal equilibrium. The sublimation rate is expressed in meters per second, representing 
the rate at which an icy surface would retreat under these conditions. We assume two 
cases: isothermal grains where the heat absorbed is radiated isotropically back into space 
("cold grains"), and grains at the temperature of the sub-solar point ("hot grains"). The 
sublimation rate is strongly dependent on the Bond albedo, which is the fraction of total 
power incident on a body that is reflected from the surface. High-albedo pure ice grains have 
low sublimation rates since they absorb little heat and remain cold. However, adding even 
a small amount of refractory material to the ice lowers the albedo significantly, resulting in 



greater absorption of heat and an increase in the temperature of the grain (Clark 1982). 



The radius of the ice grain halo, L, is related to the grain velocity, v{a) (m s ^), radius, 
a (m), and sublimation rate, a (m s^^), by: 
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via) 



(6) 



We utilized the velocity-radius relationship given in Equation ^ and measured L by 
comparing the surface brightness profile on each night to the model profile. We defined the 
edge of the halo as where the brightness fell to within 10% of the model. We acknowledge 
that this approach is model-dependent but use it to explore the concept of an ice grain halo 
in this context. The radius determined varied between ~ 2900 km and 4900 km (Table |3]). 



The median radius was 3450 km, consistent with results from Yang et al. (2009) who 



found that ice grains persisted up to 3500 km from the nucleus in 2007 Oct. We plot the 
radius-albedo relations for hot and cold ice grains that can survive across a 3450 km radius 
aperture in Figure Cold ice grains with a radius of 1 /im (c.f. Yang et al. 2009) must 
have a Bond albedo of ~ 0.35 in order to survive, while micron-sized grains in the high 
temperature limit must have a higher Bond albedo, of ~ 0.85. Of these two approximations, 
the cold grain case is probably more relevant to ice grains in 17P/Holmes, given their small 
size and the likelihood that they were rotating rapidly. 



The results of Ishiguro et al. (2010) suggest that the geometric albedo of the coma was 



between ~ 0.05 and 0.20. We note that the results from Ishiguro et al. (2010) are based on 



measurement of the flux reflected from all material in a 5300 km radius aperture, which 
would have included a significant contribution from lower-albedo dust grains as well as 
high-albedo ice grains. In this sense, our result that the Bond albedo of ice grains in the 
halo must have been ~ 0.35 appears to be consistent with their findings. 

Comparatively little contaminating material is required to significantly darken icy 



material. Clark (1982) found that contaminating 1 mm snow grains with 0.1 yUm soot 
particles that amounted to only 10~^ of the sample weight lowered the albedo from 0.94 to 



0.40. Results from Yang et al. (2009) suggested that only ice grains with < 10% impurities 
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could survive long enough to travel 3500 km from the nucleus before sublimating. |Beer et 



al. (2006) found that micron-sized pure ice grains beyond ~ 2 AU should survive for weeks 
or even months, while dirty ice grains (10% dirt by mass) should survive for only hours, 
which agrees well with our results. 

The ice grain halo radius varied during our observing campaign, although it is difficult 



to discern any obvious trends (Figure 10). We note that the halo was larger after periods 
of known activity, such as on 2007 Nov. 6 - our observation closest to the original outburst 
- and on 2007 Nov. 12 after the much smaller second outburst. The radius appeared to 
decrease after these periods of activity until the heliocentric distance of 17P/Holmes began 
to change appreciably. As the heliocentric distance increased from 2.63 AU on 2007 Dec. 11 
to 2.90 AU on 2008 Feb. 12, the temperature of the ice grains must have dropped, causing 
a decrease in sublimation rates and thus, an increase in the survival time of the grains. 
We attempted to fit the increase of the ice grain halo radius using the thermal equilibrium 



model described previously. Figure ( 10 ) shows the data with results from models of the halo 
radii expected during our observations, assuming the halo consists of isothermal grains. 
The increase in heliocentric distance by itself cannot account for the halo size variations. 
We suggest that changes in the level of activity, size and albedo of the ice grains ejected, or 
the grain velocities may also play important roles in the extent of the ice grain halo. 

Finally, we compare the expected mass loss rate from a sublimating ice grain halo to 
the observed production rates of OH molecules using aperture photometry from this work 



and production rates taken from Schleicher (2009). The production rates were derived from 



narrowband photometry of 17P/Holmes' inner coma between UT 2007 Nov. 1 and 2008 
Mar. 5. Since the published production rates were corrected for 10^ km radius apertures, 
we also use 10^ km radius apertures. We convert the apparent magnitudes to geometric 



scattering cross-sections using Equation ffl and assumptions detailed in Section 4.1, To 
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estimate the mass loss rate expected from a sublimating ice grain halo, we multiply the 
cross-section of material (presumed to be 100% ice) by the density of the ice (assumed 
to be 1000 kg m~^) and the sublimation rate (m s~^) for ice grains at the corresponding 



heliocentric distances of 17P/Holmes. Figure 11 compares the expected mass loss rates for 



ice grain halos composed of grains with albedos of and 0.3 to the derived production rates 



from Schleicher (2009). We find that the mass loss rate for an ice halo composed of grains 
with albedo = 0.3 drops off more quickly than the derived water production rate. The 
changing mass loss rate for an ice grain halo with grains of albedo = matches the trend 
in the water production rate well but the albedo conflicts with our previous determination 
that the extent of the halo is best modeled by sublimating ice grains with albedos of ~ 0.2 
- 0.4. We acknowledge that the model is underconstrained given that the ice grains may 
not be isothermal, would have a range of sizes and albedos, and those that determine the 
observed radius of the ice grain halo may not be the grains that supply the sublimated gas. 
We conclude by simply saying that the decrease in water production rates agrees loosely 
with the expected decrease in mass loss rates for an ice grain halo consisting of low-albedo 
ice grains. 



5. Summary 

We have used aperture photometry to examine the inner coma of 17P/Holmes between 
UT 2007 Nov. 6 and 2008 Feb. 12 following its massive outburst in 2007 Oct. Our results 
include: 

1. The nucleus of 17P/Holmes remained active for at least several months following its 
spectacular outburst in 2007 Oct., although at a level very low compared to the peak 
outburst rates. The inner coma decreased in brightness by a factor of ~ 50 during our 
observations, but an examination of aperture-crossing times shows that fresh material 
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must have been continually supplied from active regions on the nucleus. 

2. A second much smaller outburst occurred on UT 2007 Nov. 12, ejecting approximately 
8 X 10^ kg of material that had traveled up to 5000 km by the time of our observations. 
We propose that a fragment broke away from the surface and rapidly broke apart, 
releasing a significant amount of dust into the coma. 

3. The surface brightness profile of 17P/Holmes had a measured slope between -0.2 and 
-0.3. This is less steep than expected for a steady-state coma consistent with the 
fragmentation of icy grain aggregates. 

4. A halo consisting of freshly released ice grains persisted around the nucleus, extending 
to radii of several thousand km. This halo was likely responsible for the anomalously 



high water production rates observed by Schleicher (2009) several months after the 



outburst, although the bulk of gas and dust production could originate from distinctly 
different populations of grains. 

5. The radius of the ice grain halo is compatible with an origin in isothermal micron-sized 
grains having a Bond albedo of ~ 0.35. The ice grains may have been contaminated 
with small amounts of refractory material. 



Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT 
and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated 
by the National Research Council (NRC) of Canada, the Institut National des Science 
de rUnivers of the Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products produced at the 
TERAPIX data center located at the Institut d'Astrophysique de Paris. We appreciate 
grant support from the NASA Outer Planets Research Program to DJ. RS acknowledges 
support from the NASA Postdoctoral Fellowship Program. Part of the research was carried 



- 20 - 



out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract 
with the National Aeronautics and Space Administration. 
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Table 1. Journal of Observations 



TIT Datp 


T^'iVnri^^iTrp f imp \^\ 




A fAUl 2 


rv \f\ PO"! ^ 


Spfllp rinn/niYpll 

kJ^CllC liVlll/ ^lA-ClJ 


2007 Nov 6.4 


5 


2.49 


1.62 


13.7 


218 


2007 Nov 8.5 


5 


2.50 


1.62 


13.3 


218 


2007 Nov 9.5 


5 


2.50 


1.62 


13.1 


218 


2007 Nov 10.5 


5 


2.50 


1.62 


12.9 


218 


2007 Nov 11.5 


50 


2.51 


1.62 


12.7 


218 


2007 Nov 12.5 


5 


2.51 


1.62 


12.5 


219 


2007 Nov 13.5 


50 


2.52 


1.63 


12.3 


219 


2007 Nov 14.5 


50 


2.52 


1.63 


12.2 


219 


2007 Nov 15.5 


50 


2.52 


1.63 


12.0 


219 


2007 Dec 11.5 


70 


2.63 


1.77 


12.5 


238 


2007 Dec 13.4 


140 


2.63 


1.78 


12.8 


239 


2008 Jan 3.4 


200 


2.73 


2.01 


16.5 


270 


2008 Jan 12.3 


200 


2.77 


2.13 


17.8 


287 


2008 Jan 16.4 


200 


2.78 


2.19 


18.3 


295 


2008 Feb 12.3 


50 


2.90 


2.62 


19.8 


352 



Heliocentric distance 
^Geocentric distance 
^Phase angle 
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Table 2. Characteristics of Coma Material 



UT Date 


m ^ 






Mass 


Mass Loss Rate 








[10^ m2] 


[10^ kg] 


[kg s-i] 


2007 Nov. 6.4 


11.7 


8.1 


166.0 


221.4 


4870 


2007 Nov. 8.5 


12.3 


8.7 


95.9 


127.8 


2810 


2007 Nov. 9.5 


12.5 


8.9 


79.7 


106.3 


2340 


2007 Nov. 10.5 


12.7 


9.1 


68.2 


90.9 


2000 


2007 Nov. 11.5 


12.8 


9.2 


59.1 


78.8 


1730 


2007 Nov. 12.5 


12.6 


9.0 


71.1 


94.8 


2090 


2007 Nov. 13.5 


13.0 


9.4 


51.5 


68.6 


1510 


2007 Nov. 14.5 


13.2 


9.6 


43.3 


57.7 


1270 


2007 Nov. 15.5 


13.3 


9.7 


39.3 


52.4 


1150 


2007 Dec. 11.5 


15.0 


11.1 


10.4 


13.9 


310 


2007 Dec. 13.4 


15.2 


11.2 


9.3 


12.4 


270 


2008 Jan. 3.4 


16.1 


11.7 


6.4 


8.5 


190 


2008 Jan. 12.3 


16.5 


11.9 


5.1 


6.7 


150 


2008 Jan. 16.4 


16.6 


11.9 


5.2 


6.9 


150 


2008 Feb. 12.3 


17.6 


12.4 


3.1 


4.2 


90 



Apparent magnitude within 2500 km projected radius aperture 
^Absolute magnitude within 2500 km projected radius aperture 
^Cross-section of material within 2500 km projected radius aperture 
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Table 3. Surface Brightness Profile Characteristics 



Date (UT) Surface Brightness Profile Gradient Ice Grain Halo Radius [km] 



2007 Nov. 6.4 


-0.25 


4870 


2007 Nov 8 5 


-0.26 


3680 


2007 Nov. 9.5 


-0.26 


3310 


2007 Nov. 10.5 


-0.25 


3120 


2007 Nov. 11.5 


-0.25 


3080 


2007 Nov. 12.5 


-0.25 


4590 


2007 Nov. 13.5 


-0.28 


3410 


2007 Nov. 14.5 


-0.26 


3300 


2007 Nov. 15.5 


-0.26 


2940 


2007 Dec. 11.5 


-0.26 


3250 


2007 Dec. 13.4 


-0.25 


3730 


2008 Jan. 3.4 


-0.25 


4160 


2008 Jan. 12.3 


-0.25 


4240 


2008 Jan. 16.4 


-0.24 


4330 


2008 Feb. 12.3 


-0.18 


4500 
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2007 Nov. 6 1 2007 Nov. 8 I 2007 Nov. 9 1 2007 Nov. 10 




2008 Jan. 12. 1 2008 Jan. 16 I 2008 Feb. 12 




10^ km ^ I lO'km -J- I 10° km 



Fig. 1. — The full MegaCam field of view is shown for each observation of 17P/Holmes. 
Each image is oriented with north up, east left. The anti-solar vector, as projected on the 
plane of the sky, is denoted by a white arrow. The final image is slightly larger due to the 
use of a larger dithering patten. As the coma surrounding the nucleus expands the surface 
brightness drops and the coma becomes more diffuse. 
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Fig. 2. — The inner coma of 17P/Holmes as seen throughout our observations. Each image 
is 1' X 1' and is oriented with north up and east left. The anti-solar vector, as projected on 
the plane of the sky, is denoted by a white arrow. 
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2007 Nov. 6 


2007 Nov. 8 


2007 Nov. 9 


2007 Nov. 10 





















2007 Nov. 11 I 2007 Nov. 12 I 2007 Nov. 13 ■ 2007 Nov. 14 




2007 Nov. 15 12007 060.11 12007 060.13 12008 0311.3 




2008 Jan. 1 2 ■ 2008 Jan. 1 6 ■ 2008 F6b. 1 2 




Fig. 3. — The central arcminute region of each image centered on the nucleus in north- 
up, east-left orientation with the projected anti-solar vector marked by a white arrow. The 
isophotes are of arbitrary brightness and were chosen to highlight the morphological evolution 
of the coma. The innermost coma was initially elongated but became more symmetrical 
towards the end of our observations. 
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2007/1 1 /06 2007/1 2/26 2008/02/1 4 

Date of Observation (UT} 



Fig. 4. — The fading of 17P/Holmes over four months as observed using an aperture with a 
projected radius of 2500 km. The absolute magnitudes were calculated using Equation ([T]) 
and have been fitted with an exponential function. The brightness of the inner coma dropped 
by a factor of ~ 50 during our observing campaign. 
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10.0L_^ , , : I 1^ , , ^ J 

2007/1 1 /06 2007/1 1 /1 1 2007/1 1 /1 6 
Date of Observation (UT) 

Fig. 5. — The photometric evolution of the inner coma of 17P/Holmes during 2007 Nov. as 
measured with a 2500 km radius aperture. The best fit to the data was found by fitting 
an exponential function and is shown by the solid line. While there was an overall decrease 
in the absolute magnitude of 17P/Holmes during this time, an apparent brightening can be 
seen to occur on 2007 Nov. 12. The deviation of the observed magnitude from the value 
expected according to the best fit is 0.3 mag, corresponding to a 3-cr deviation. 
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1000 10000 
Projected Aperture Radius [km] 



Fig. 6. — The surface brightness profiles of the inner coma, centered on the nucleus, from 
2007 Nov. 6 UT (top) to 2008 Feb. 12 UT (bottom). Error bars are plotted and represent a 
10% uncertainty in the sky background estimate. Each profile has a visible bump near the 
center and has a slope of ~ -0.2 to -0.3 between the projected distances of 10000 km and 
25000 km. The surface brightness of the inner coma, as observed at 25000 km, decreases 
monotonically with time. The surface brightness profile within ~ 5000 km on 2007 Nov. 
12 is brighter than expected, which is consistent with a second outburst occurring at the 
nucleus shortly before our observation. 



33 




Fig. 7. — The estimated cross-section and mass of material within 2500 km of the nucleus, 
calculated using Equation Q and the assumption that the material consists of dust grains 
with radii of 1 /im with an albedo of 0.1 and a density of 1000 kg m~^. 
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1000 10000 1000 10000 1000 10000 

Projected Aperture Radius [km] 



Fig. 8. — The surface brightness profile of 17P/Holmes' inner coma (sohd hne) compared 
to a generated profile (dashed hne). The generated profile consists of an unresolved central 
source with coma that has the same slope as the 17P/Holmes coma between 10000 km and 
25000 km. Both profiles have been convolved with a 1" Gaussian function to eliminate the 
effects of variable seeing from night to night. The generated profile fails to simulate enough 
flux in the innermost coma of 17P/ Holmes, leading to the possibility that an ice grain halo 
persisted at small nucleo-centric distances. 
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-|q6| , , , 1 , , r 1 , , , 1 , , , 1 , , r 




Albedo 

Fig. 9. — Grain radius versus albedo relation for ice grains having a sublimation lifetime 
sufficient to travel 3450 km - the median of our estimated ice grain halo radii. The dashed 
and solid lines depict results for low and high grain temperature models, as described in 
the text. Cold micron-sized ice grains must have had an albedo of ~ 0.35 (dotted lines) 
in order to move 3450 km from the nucleus at 550 m s~^, which suggests that the ice was 
contaminated with albedo-lowering regolith and was not pure. 
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2007/11/01 2007/11/28 2007/12/25 2008/01/21 2008/02/17 

Date of Observation [UT] 

Fig. 10. — The radius of the observed ice grain halo (filled circles, from Figure 8) com- 
pared to the expected radius of an ice grain halo given the varying heliocentric distance of 
17P/Holmes. At larger heliocentric distances, and correspondingly lower temperatures, ice 
grains can survive for longer periods of time so we expect the radius of the halo to increase 
as 17P/Holmes moves away from the Sun. We show the radius expected for halos with ice 
grains of varying sizes and constant albedos of 0.3. None of our scenarios provide a good fit 
to the data. We suggest that changes in other parameters, such as the ejection velocity of 
the ice grains, or the albedo, may play an important part in the maximum extent of the ice 
grain halo. 
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Fig. 11. — A comparison of the derived water production rates in the inner 10^ km of the 
coma (crosses; 



Schleicher 2009 ) and the expected mass loss rates for a sublimating ice grain 



halo composed of grains with albedos of (filled circles) and 0.3 (triangles). The calculated 
mass loss rates have been plotted such that the mass loss rate for albedo = on UT 2007 
Nov. 11.5 coincides vertically with the data point representing the derived water production 
rate on UT 2007 Nov. 11.2. The decrease in the mass loss rate for an ice grain halo with an 
albedo of appears consistent with the decrease in water production rates, while the mass 
loss rate for a halo with higher-albedo grains would decrease more steeply. 



